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Parametric instabilities have long been studied as a potentially limiting effect in high-power 
interferometric gravitational wave detectors. Until now, however, these instabilities have never 
been observed in a kilometer-scale interferometer. In this work we describe the first observation of 
parametric instability in an Advanced LIGO detector, and the means by which it has been removed 
as a barrier to progress. 


INTRODUCTION 

Opto-mechanical interactions, the moving of masses 
by radiation pressure, are typically of such a tenuous 
nature that even carefully designed experiments may 
fail to observe them. In the extreme environment of a 
high-power interferometric gravitational wave detector, 
however, these effects arise spontaneously. This is true 
despite the fact that these instruments feature multi¬ 
kilogram optics [T], unlike the micro-gram or nano-gram 
optics generally used in experiments which target such 
effects [2H5]. 


Since the seminal work of Braginsky, Strigin, and Vy¬ 
atchanin in 2001 [6], optical parametric instabilities have 
been extensively studied as a potential limit to high- 
power operation of interferometric gravitational wave de¬ 
tectors I3S5]. These instabilities are of acute interest 
because, if left unchecked, they will limit detector sensi¬ 
tivity by limiting circulating power. 

We report on the first observation of a self-sustaining 
parametric instability in a gravitational wave detec¬ 
tor and the subsequent quenching of this instability. 
This observation, at the LIGO Livingston Observatory 
(LLO) [36], the culmination of more than a decade of the- 
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FIG. 1. Parametric instabilities transfer energy from the arm 
cavity field to a mechanical mode of an interferometer test 
mass. Since the rate of energy transfer depends on the ampli¬ 
tude of the excited mechanical mode, this can be a runaway 
positive feedback process. 


oretical calculation, numerical modeling, and lab-scale 
experimentation, and serves as confirmation that models 
which have been built to understand the phenomenon of 
parametric instability are substantially correct. 

A variety of techniques have been suggested for over¬ 
coming parametric instabilities in gravitational wave de¬ 
tectors [ 37 H 43 ] . There are a few notable categories of 
mitigation techniques: avoid instability by changing the 
radius of curvature of one or more optics [371 [43], actively 
damp mechanical modes as they become excited m, and 
prevent instability by increasing the loss of the test mass 
mechanical modes [29l [39l [40l [44] . The first of these tech¬ 
niques has been demonstrated and found to be effective 
at LLO. 


THEORY OF PARAMETRIC INSTABILITY 


Parametric instabilities (PI) operate by transferring 
energy from the fundamental optical mode of the inter¬ 
ferometer, which with nearly 1 MW of circulating power 
can be as much as 40 J, into an interferometer optic’s 
mechanical mode (see figure]^. Energy transfer takes 
place via the radiation pressure driven opto-mechanical 
interaction and the modulation of the fundamental field 
by the excited mechanical mode [25] (see figure . 

The parametric gain for a given test mass mechanical 
mode is given by 


^ _ STrQm^arm ^ ^2 
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(using the notation of reference [25], and correcting an 
error therein [45]). The fixed parameters used in equation 


[2 are: c is the speed of light, A the laser wavelength, M 
the mass of the optic, uJrn the angular resonant frequency 
of mechanical mode with index m, and Bm,n the overlap 
between mechanical mode m and optical mode n. 

Clearly, the potential for instability grows with in¬ 
creasing power, since the parametric gain of all modes 
is proportional to the power circulating in the interfer¬ 
ometer arm cavities, Parm- Active and passive damping 
techniques for defusing PI operate by reducing the Q of 
an otherwise unstable mechanical mode, while avoidance 
techniques work by changing the optical gain Gn to pre¬ 
vent instability. 

The e-fold growth time for an unstable mode is given 
by 

2Qm 

- 1 ) 

Note that for Rm = 0 this gives a negative value equal 
to the usual decay time of a mechanical mode with fre¬ 
quency fm = ujrnl‘^7r and quality factor Qm- For values 
of Rjn > 1 equation [^ gives a positive value, indicating 
exponential growth. 

The circulating power level Parm at which a mode be¬ 
comes unstable is referred to as the threshold power Pm 
for that mode, and is found by rearranging equation [^ 
with Rm = 1 and Parm = Pm- At threshold the opto¬ 
mechanical interaction is putting energy into the mode 
at the same rate that it is being dissipated, so ^ oc. 

OBSERVED PARAMETRIC INSTABILITY 

Parametric instability was first observed in the Ad¬ 
vanced LIGO detector recently installed at LLO (see ref¬ 
erence [46] for detailed information about the detectors). 
The instability grew until it polluted the primary gravi¬ 
tational wave output of the detector by aliasing into the 
detection band and saturating detection electronics (see 
figure [^. The time scale for growth was long enough 
to allow for manual intervention (a reduction of the laser 
power input to the interferometer) before control was lost 
due to this saturation of the readout electronics. 

The test mass (TM) mechanical mode responsible 
for the instability was identified as the 15.54 kHz mode 
shown in figure [^ The higher-order mode spacing of the 
Advanced LIGO arm cavities is 5.1 ± 0.3 kHz, and the 
optical resonance width is 80 Hz, such that a 3^^ order 
transverse optical mode can provide energy transfer from 
the fundamental optical mode to this mechanical mode 
(see figure and 1^ . 

By measuring the e-folding growth and decay time of 
the excited acoustic mode as a function of circulating 
power in the interferometer arm cavities, we can com¬ 
pute the parametric gain R and mechanical mode qual¬ 
ity factor Q as given in equation [^ Our measurements 
were performed by operating the interferometer above 
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FIG. 2. The excitation due to parametric instability was first 
observed as an exponentially growing feature in the detector’s 
primary output (the gravitational wave channel). The feature 
is aliased into the detector band by the 16 384 Hz sample rate 
of the digital system which causes it to appear at 845 Hz. 
The growing instability eventually causes saturation of the 
electronic readout chain, which appears as broadband con¬ 
tamination of the detector output channel (visible between 
2500 and 3000 s). In the above data, the power into the in¬ 
terferometer was decreased before saturation caused the in¬ 
terferometer control systems to fail (see figure |^. 


FIG. 3. The test mass mechanical mode and cavity optical 
mode responsible for the observed parametric instability are 
shown. The left panel is the test mass front surface displace¬ 
ment due to the mechanical mode, while the right panel is the 
radiation pressure induced displacement (red is positive and 
blue negative in both panels). Both of these modes occur at 
15.5 kHz and they have an overlap factor of Bm,n = 0.1. 

the threshold power to excite the TM mechanical mode 
and then reducing the power below the instability thresh¬ 
old and watching the mode amplitude decrease, as shown 
in figure 

The observed unstable mode has R = 2 with Parm = 
50 kW, and the associated mechanical mode has Q = 
12 X 10^. This Q-factor is in the range expected given 
similar measurements of Advanced LIGO test masses [42l 
[47], and the parametric gain is as predicted by equation 
[^for a high, but far from maximal value of Gn- While 
the 90% confidence limit for this mode is P = 11, as seen 
in figure 5% of simulated values are higher than the 
observed value, and 1% have R > 100. 
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FIG. 4. The amplitude of the excitation shown in figure is 
fitted at times with different values of Parm to find its growth 
(or decay) time-scale as a function of power. The growth 
of the PI is clearly visible above the noise after 1000 s, with 
an e-folding growth time of 240 s, until 2000 s at which point 
the readout electronics begin to saturate. A little more than 
4000 s into the plotted data, the excited mechanical mode 
is seen decaying with r = —900 s. According to equation 

these data imply a threshold power of 25 kW and Q = 
12 X 10® for this mode. 

DEFUSING PARAMETRIC INSTABILITY 

Parametric instability depends on several potentially 
modifiable features of a gravitational wave detector, two 
of which can be taken advantage of in Advanced LIGO 
without modifications to the interferometer core optics. 
First, instability requires coincident resonant frequencies 
of a test mass mechanical mode and an arm cavity op¬ 
tical mode (i.e., Gn in equation must be large at the 
mechanical mode frequency ujm for an optical mode with 
non-vanishing overlap 

The coincidence of resonant frequencies can be modi¬ 
fied by changing the radius of curvature (RoC) of one of 
the test masses in the cavity affected by PI, through its 
effect on the transverse mode spacing. Advanced LIGO 
arm cavities use mirrors with ^ 2 km RoC, so for the 
3^^ order optical mode shown in figure for instance, 
a change of 1 m is sufficient to change its resonant fre¬ 
quency by 80 Hz, or one cavity line width. 

A second approach to defusing PI comes from the weak 
radiation pressure coupling of energy from the fundamen¬ 
tal optical mode to the mechanical modes of the test 
mass, which must be sufficient to overcome the energy 
loss from the mechanical mode. This can be seen in the 
linear dependence of Rm on Qm in equation The qual¬ 
ity factor of mechanical modes in the Advanced LIGO 
test mass optics is roughly 10^, and as such is easily 
spoiled. Since the highest parametric gain likely to be 
seen in Advanced LIGO is ^ 10, reducing Qm to less than 
10^ for potentially unstable modes would suffice (see fig- 
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FIG. 5. The maximum gain expected at full power, Parm = 
800 kW, for all potentially unstable modes (90% confidence 
value for a single test mass). Since the mechanical mode fre¬ 
quencies and optics parameters are known with limited preci¬ 
sion, the parametric gain of a particular mode cannot be com¬ 
puted exactly, and Monte Carlo methods must be used [25]. 
The “90% confidence maximum” shown here is the value for 
which 90% of the 342,881 Monte Carlo computations gave a 
lower gain. Notably, the observed unstable mode is the mode 
with the highest predicted parametric gain (marked with a 
green star), and the observed gain is about a factor of 3 higher 
than the value shown here (P = 2 at 50 kW, or P = 32 at 
800 kW). This is not so much bad luck as a trials factor; there 
are 4 test masses, and many potentially unstable modes, so it 
is not unusual that at least one is above the 90% confidence 
limit. In fact, the statistics shown in figure [^confirm that the 
observation of at least one unstable mode was to be expected. 

ure[^. 

In anticipation of PI, all Advanced LIGO test masses 
were outfitted with electro-static actuators capable of 
damping mechanical modes associated with PI |42| . This 
method of damping TM mechanical modes was demon¬ 
strated at MIT with a prototype Advanced LIGO optic, 
and will be used to damp PI in Advanced LIGO as nec¬ 
essary. 

As discussed in the following section, the combination 
of the above techniques will likely prove sufficient if a 
manageable number of unstable modes are encountered. 
If, however, the number of unstable modes is very large, 
a broadband reduction of mechanical Q-factors may be 
required. Implementing a broadband PI control strategy, 
as described in m, will certainly prove successful, but 
it may cause some delay in the approach of Advanced 
LIGO to its design operating power. 


IMPLICATIONS FOR THE FUTURE OF 
ADVANCED LIGO 

Shortly after the first PI was observed at LLO, heating 
elements, known as “ring heaters”, were used to change 
the mirror RoC and avoid the instability. Ring heaters 


P =100kW P =800kW 

arm arm 
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FIG. 6. The number of unstable modes as a function of change 
in arm cavity mirror radius of curvature (RoG), as indicated 
by the color scale to the right of each plot. With 100 kW 
of circulating power, as expected for the first observing run, 
small changes in mirror RoG can be used to avoid parametric 
instability. Higher power levels, however, increase paramet¬ 
ric gain; note the different color scales used on the left and 
right panels. At Advanced LIGO’s design operating power 
of Parm = 800 kW RoG adjustments can at best be used to 
reduce the number of unstable modes. 


were included in the Advanced LIGO design to compen¬ 
sate for RoC changes due to absorption of the fundamen¬ 
tal optical mode, and can change the RoC of the optics by 
several 10s of meters. After adjusting the ring heaters to 
produce roughly 2 m of RoC change of the arm-cavity end 
mirrors, the parametric gain of the observed instability 
at 15.53 kHz was reduced below unity and the interfer¬ 
ometer was operated with nearly 100 kW of circulating 
power for more than 12 hours. 

Despite this success, it must be recognized that many 
mechanical modes of the interferometer test masses can 
be driven to instability; figure shows the maximum ex¬ 
pected parametric gain of all potentially unstable modes. 
From the left panel of figure we can conclude that us¬ 
ing the ring heaters to find a Pl-free zone is likely to 
succeed at present (with Parm ^ 100 kW), but the right 
panel tells us that at the design value of Parm = 800 kW 
RoC changes will not be sufficient. 

Several major periods of data taking are expected 
with Advanced LIGO at power levels below the design 
power |48|. Figure]^ shows the number of modes that are 
likely to be unstable as a function of circulating power 
in the interferometer. At the final operating power, for 
instance, more than 40 modes are likely to be unstable. 

In order to operate at the design power level Advanced 
LIGO will likely need to employ multiple methods of de¬ 
fusing PL changing the RoC to reduce the number of un¬ 
stable modes; actively damping the remaining unstable 
modes with electrostatic actuators; and possibly reducing 
the Q of the test mass modes with passive dampers. 
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FIG. 7. The probabilities of having various numbers of un¬ 
stable modes anywhere in the interferometer (with 4 test 
masses), as a function of circulating power, Parm- Each curve 
gives the probability that at least N = {1,2, 3, 5,...} unsta¬ 
ble modes are observed in an Advanced LIGO interferome¬ 
ter at a given level of circulating power in the arm cavities. 
The vertical grey lines mark the threshold power of the mode 
first observed to be unstable, near 25 kW of circulating power, 
the power level at which the first observing run of Advanced 
LIGO is planned, roughly 100 kW, and the design power level 
of 800 kW. These data indicate that more than 40 modes will 
likely need to be damped or otherwise defused in order to 
operate at 800 kW. 

CONCLUSIONS 

Parametric instabilities, studied in great depth and 
feared as a limitation to the attainable operating pow¬ 
ers of interferometric gravitational wave detectors, have 
been observed for the first time in an Advanced LIGO 
detector. The behavior of the observed instability was 
found to be largely in agreement with models of the ef¬ 
fect, implying that no significant ingredients have been 
omitted in the theoretical analysis. 

Furthermore, the observed PI has been quenched by 
thermally tuning the optical resonance of the interferom¬ 
eter away from the resonance of the associated mechan¬ 
ical mode. This approach to PI, while sufficient for a 
small number of potentially unstable modes, may not be 
sufficient at higher operating powers where many modes 
are available for runaway excitation. 

Thanks to many years of theoretical and experimental 
work on parametric instabilities, now informed by the ob¬ 
servations described in this paper, the challenge faced by 
high-power interferometric gravitational wave detectors 
is clear and well understood. While the necessary miti¬ 
gation techniques are not trivial, a suitable combination 
of thermal tuning, active damping of excited mechan¬ 
ical modes, and passive reduction of mechanical mode 
Q-factors is expected to be sufficient to allow Advanced 
LIGO to operate stably at full power. 


The authors would like to acknowledge the extensive 
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Moscow State University colleagues Vladimir Braginsky, 
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